KCrF 3 : Electronic Structure, Magnetic and Orbital Ordering from First Principles 



o 
o 

(N 



^1- 



a 

C 

o 
o 



> 
in 

in 
cn 

i> 
o 

o 



Gianluca Giovannetti 1,2 , Serena Margadonna 3 , Jeroen van den Brink 1 ' 4 
1 Institute Lorentz for Theoretical Physics, Leiden University, 
P.O. Box 9506, 2300 RA Leiden, The Netherlands 
2 Faculty of Science and Technology and MESA+ Research Institute, 
University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands 
3 School of Chemistry, University of Edinburgh, West Mains Road, Edinburgh EH9 3JJ, UK 
Institute for Molecules and Materials, Radboud Universiteit Nijmegen, 
P.O. Box 9010, 6500 GL Nijmegen, The Netherlands. 
(Dated: February 1, 2008) 

The electronic, magnetic and orbital structures of KCrF3 are determined in all its recently iden- 
tified crystallographic phases (cubic, tetragonal, and monoclinic) with a set of ab initio LSDA 
and LSDA+U calculations. The high-temperature undistorted cubic phase is metallic within the 
LSDA, but at the LSDA+U level it is a Mott insulator with a gap of 1.72 eV. The tetragonal and 
monoclinic phases of KCrF3 exhibit cooperative Jahn- Teller distortions concomitant with staggered 
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orbital order. We find that the energy gain due to the Jahn- Teller distortion is 



82/104 meV per chromium ion in the tetragonal/monoclinic phase, respectively. These phases show 
A-type magnetic ordering and have a bandgap of 2.48 eV. In this Mott insulating state KCrF3 has 
a substantial conduction bandwidth of 2.1 eV, leading to the possibility for the kinetic energy of 
charge carriers in electron- or hole-doped derivatives of KCrF3 to overcome the polaron localization 
at low temperatures, in analogy with the situation encountered in the colossal magnetoresistive 
manganites. 



PACS numbers: 71.45.Gm, 71.10.Ca, 71.10.-w, 73.21.-b 



INTRODUCTION 

About a decade ago, the discovery of the Colossal Mag- 
netoresistance (CMR) effect in doped manganites caused 
a surge of interest in these perovskite oxides 0, Q ■ The 
particular physical properties of the CMR materials are 
related to the fact that their parent compound LaMn03 
contains Mn 3+ ions with 4 electrons in its d-shell. On the 
one hand, the presence of these Jahn- Teller active ions 
leads to a strong coupling between the electrons and the 
lattice, giving rise to polaron formation which is widely 
perceived to be essential for the CMR effect [H, 0| • On 
the other hand, when doped, the d A high spin state leads 
via the double exchange mechanism to a ferromagnetic 
metallic state with a large magnetic moment making the 
system easily susceptible to externally applied magnetic 
fields Moreover the presence of strong electron cor- 
relations and an orbital degree of freedom, to which the 
Jahn- Teller effect is directly related, leads to an extraor- 
dinary rich phase diagram at higher doping concentra- 
tion, displaying a wealth of spin, charge, orbital and mag- 
netically ordered phases It follows that the high 
spin d A state of Mn 3+ is intimately related to a plethora 
of physical phases, effects and properties. However, it 
is important to note that the high spin <i 4 state is not 
exclusive to trivalent manganese. 

Formally high spin Cr 2+ is electronically equivalent 
to Mn 3+ . However, divalent chromium due to its low 
ionization potential is rarely found in solid state sys- 
tems, making KCrF3 an intriguing example. Recently, we 
characterized in detail the temperature-dependent crys- 
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FIG. 1: (a) Cubic perovskite crystal structure of KCrF3 at 
high temperature, (b) Schematic representation of the inter- 
mediate temperature tetragonal structure with CrF6 octahe- 
dra elongated in an alternate fashion along the (110) (blue) 
and (1-10) (red) directions, (c) Schematic representation of 
the low temperature monoclinic structure with CrF6 octahe- 
dra elongated along the (1-10) (blue) and (001) (red) direc- 
tions 



tallographic phase diagram of KCrF3 revealing strong 
structural, electronic, and magnetic similarities with 
LaMn03 [H, 0], including the presence of Jahn- Teller 
distortions and orbital ordering, and orbital melting at 
high temperature. Here we report ab initio electronic 
structure calculations on the different phases of this com- 
pound, both within the local spin density approximation 
(LSDA) and the LSDA+U, in which local electron corre- 
lation effects are partially accounted for. The results of 
our calculations clearly show that KCrF 3 and LaMn03 
are not only structurally, but also electronically very sim- 
ilar. 

KCrF3 displays three different crystallographic struc- 
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tures, see Fig. [TJ At very high temperatures, the sys- 
tem is a cubic perovskite (space group Pm-3m) Be- 
low 973 K, the JT-active high spin Cr 2+ ion induces a 
lattice distortion to a body-centered tetragonal unit cell 
(space group I4/mcm), isostructural to the Cu 2+ ana- 
logue, KCuF 3 . In the tetragonal phase, the CrF 6 octa- 
hedra are distorted, leading to short Cr-F bonds along 
the c axis and alternating long-short Cr-F bonds in the 
ab plane, indicative of the presence of a staggered type 
of orbital ordering. On cooling below the 250 K, there 
is a phase transition to a more complicated monoclinic 
structure with space group I2/m. 

From our density functional calculations we find that 
the tetragonal phase of KCrF3 is a strongly correlated 
insulator with a gap of 0.49 eV in LSDA and 2.48 eV in 
LSDA+U (with U = 6 eV). For this value of U, the cal- 
culated relaxed lattice structure is in excellent agreement 
with the experimental one. In the tetragonal orbitally- 
ordered phase we find a crystal-field splitting between 
the Cr tig and e g states is 1.0 eV and a total energy 
gain related to the Jahn- Teller splitting of the e g states 
of 0.328 eV per formula unit containing four chromium 
ions. The cooperative Jahn- Teller distortion is accom- 
panied by A-type antiferromagnetic spin ordering in a 
similar fashion to LaMn0 3 [lJJ EH El, El- We find a 
magnetic moment of 3.85 ub per Cr 2+ ion, in excellent 
agreement with experiment 14| and in-plane and out- 
of-plane superexchange parameters of 2.6 meV and -3.4 
meV, respectively. 

The system displays antiferrodistortive ordering of the 
3rf3 X 2 _ r 2 and id^, y -2_ r 2 orbitals in the ab plane -an or- 
dering motif very different from the orbital ordering in 
KCUF3 that gives rise to a quasi one-dimensional spin 
chain formation, and rather resemblin g th e orbital order- 
ing in LaMnOa [IE El El EE EE E3, M El • Along the 
c axis, the orbital ordering pattern in KCrF3 is rotated 
by 90 degrees in consecutive layers. This in contrast to 
the manganite where the ordering along the c axis is a 
uniform repetition of the in-plane orbital structure. An- 
other difference with LaMnOs is that the e g bandwidth 
in the chromium compound, as computed within LSDA, 
is about a factor of two smaller. However, this is com- 
pensated in LSDA+U, which shows a bandwidth of the 
lower Hubbard e g band of 2.1 eV. 

On cooling below the 250 K, KCrF3 shows a phase 
transition to a more complicated monoclinic structure 
(space group I2/m) with four chromium atoms in the unit 
cell. Our calculations show that in the monoclinic phase 
an A-type magnetic structure is also realized and that the 
Jahn- Teller energy is lowered, leading to an even stronger 
orbital ordering. However, the resulting electronic gap 
and magnetic moment of the compound barely change. 

In the following we will present the electronic struc- 
ture calculations for the three different crystallographic 
structures. For each one, we considered several possible 
magnetic ordering structures (ferromagnetic and antifer- 
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FIG. 2: Band structure and projected density of states cal- 
culated in LSDA for the tetragonal structure of KCrF3. The 
labels on, (3i axe indicating the states 3a; 2 — r 2 , y 2 — z 2 and 
3j/ 2 — r 2 , x 2 — z 2 , respectively on neighboring Cr sites 1 and 
2 in the ab plane. 



romagnetic of A, G, and C-type) and analyzed the result- 
ing electronic properties, Jahn- Teller energies and orbital 
orderings. 



INTERMEDIATE-TEMPERATURE 
TETRAGONAL PHASE 

The structural changes which occur on lowering the 
temperature below 973 K through the cubic-to-tetragonal 
phase transition can be described in terms of two com- 
ponents: a uniform Q 3 -type tetragonal distortion, which 
shortens one lattice constant (along the c direction) and 
lengthens the other two (along the a and b directions) 
and a Q2-type staggered distortion, which introduces al- 
ternating Cr-F bond lengths in the ab plane with two 
distinct Cr-F bonds of 2.294/1.986 A. This is a textbook 
example of a cooperative Jahn- Teller distortion. The lat- 
tice parameters of the resulting body-centered tetragonal 
unit cell at room temperature are a= 6.05230 A and c— 
8.02198 A 0. 

All our self-consistent calculations are done within 
LSDA and LSDA+U using the Vienna Ab-initio Simula- 
tion Package (VASP)[I3]. Total energies for the tetrago- 
nal structure were calculated with a kinetic cutoff energy 
of 500 eV and the tetrahedron with Blochl correction, 
using 105 irreducible k-points. 



LSDA Electronic Structure of Tetragonal KCrF3 

We start our study of the tetragonal structure of 
KCrF3 at the LSDA level and then proceed to also in- 
clude local correlations within LSDA+U. We find that 
the A-type antiferromagnetic spin ordered structure is 
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FIG. 3: Contour plot of charge density corresponding to the 
occupied e g bands within the LSDA for the tetragonal struc- 
ture of KCrF3. The orbital ordering pattern is clearly seen 
along the bonds connecting the Cr ions, with 3d 3x 2_ r 2 and 
3d 3y 2_ r 2 alternating in the ab plane. The left and right panel 
are two cuts on consecutive planes along the c direction. 



the ground state. In Fig. [5J the band structure and the 
(projected) density of states are shown. The system is 
insulating with an energy gap of 0.49 eV, which is in- 
duced by the Jahn- Teller splitting of the e g states. In 
accordance with Hund's rule, the Cr 2+ ions are in a high 
spin t2 g 3 e g 1 state, giving rise to a magnetic moment of 
3.59 \ib per Cr. 

The Fermi level lies just above the bands with ti g and 
e g character, in agreement with the high spin state of the 
Cr ions. The exchange splitting is about 2.6 eV, which 
moves the minority-spin bands far above the Fermi level. 
The t2g-e g crystal field splitting, Acf is about 1.0 eV. 
The occupied Cr bands show little dispersion along the 
T—Z direction and are therefore of quasi- two-dimensional 
character, which is due to the specific ordering of the e g 
orbitals that maximizes hybridization in-plane and min- 
imizes the out-of-plane dispersion. The character of oc- 
cupied e g bands is mixed between the two types of e g 



states but is mainly coming from 3x 



bitals on neighboring Cr ions. This becomes immediately 
clear from the contour plot of charge density correspond- 
ing to the e g bands below the Fermi level, shown in Fig. [3] 
The orbital ordering is found to be staggered along the 
c direction. The unoccupied orbitals have 3d y 2_ z 2 and 
3d x 2_ z 2 character, respectively. The t2 g projected den- 
sity of states resolved for orbital character shows that 
the xy states have a different distribution in energy from 
the two-fold degenerate orbitals of yz, zx character, in 
agreement with crystal field symmetry expectations. 

The band width of the t2 g bands is 0.65 eV, while 
the Jahn- Teller split e g x bands just below and above 
the Fermi level each have a width of about 0.65 eV, al- 
most a factor of two smaller than the value of 1.2 eV in 
LaMn03 12j. However, the inclusion of local correlations 



within LSDA+U changes this bandwidth significantly. 

LSDA+U Electronic Structure of Tetragonal KCrF 3 

It is well known that the incorporation of local 
Coulomb interactions is essential to understand the phys- 
ical properties of transition metal compounds [2(. In 
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FIG. 4: Deviation of relaxed Cr-F distances dl, d2 and d3 
from their experimental as a function of the Coulomb inter- 
action U. 



LSDA+U the electron-electron interaction is dealt with 
on a meanfield level and we repeated the LSDA calcula- 
tions above within this scheme. 

We performed calculations for a series of values 
of the on-site Coulomb parameter U, namely U = 
2.0, 4.0, 6.0, 8.0 eV, adopting a value for Hund's exchange 
of Jh = 0.88 eV. In practice the exact definition of U in 
a solid is not trivial. The value that is found for this 
parameter depends on for instance the precise choice of 
the orbitals that are used in the calculation 2(J 21, 22j |. 
In order to determine its value, we performed a struc- 



tural optimization as a function of U and subsequently 
stay with the value for U for which we find an equilib- 
rium structure that matches the experimental one. This 
scheme to extract the Coulomb parameter is viable be- 
cause the on-site Hubbard U determines for a large part 
the orbital polarization of the e g states, which in turn 
causes the structural Jahn- Teller lattice distortion [l7j ]. 

Hund's exchange parameter Jh, in contrast, represents 
a local multipolc and is only very weakly screened in the 
solid and therefore close to its bare atomic value. For it 
we used the value for a high spin d 4 configuration deter- 
mined by constrained density functional calculations [13] . 
At any rate small changes of Jjj will not affect the re- 
sults of LSDA+U significantly, as U is the dominating 
parameter. 

To determine U we optimized the three inequivalent 
Cr-F distances in the tetragonal unit cell, d\, di and d^, 
while fixing the lattice parameters a, b and c by mini- 
mizing the total energy. The results are shown in Fig. [4] 
For U = 6.0 eV we find that the computed structure 
is very close to the one obtained experimentally, moti- 
vating us to adopt this value as the most reliable one. 
The obtained U value agrees well with that for Mn 4+ 
m LaMnOa (8.0 eV) calculated within a constrained 
LDA+U scheme. In KCUF3 a value of 7.5 eV is similarly 
found for Cu 2+ [ljj. As one expects that the ionic core 
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LSDA+U tetragonal structure Cr e g 



Crt 



TABLE I: Energy difference (eV) between the magnetic 
ground state (A-type) and other magnetic orderings, band 
gap A (eV) and magnetic moment (fit) of the Cr-ions in 
KCrF 3 . 
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FIG. 5: Band structure calculated in LSDA+U with U = 6.0 
and J = 0.88 eV for the tetragonal structure of KCrF3 and 
(projected) densities of states. 



potential of Cr 2+ causes the d electrons to be less local- 
ized with respect to both examples above, it is reasonable 
to find the smaller value of U — 6.0 eV for K&F3. 

Below 46 K antifcrromagnctic spin ordering is ob- 
served [8]. To study the magnetic exchange couplings 
between the Cr ions we calculated the total energy of var- 
ious magnetic structures. The different magnetic struc- 
tures we considered are: A-type (the spins are parallel 
in the ab planes and the spins are antiparallel along the 
c-axis), F-type (all spins parallel), C-type (each spin is 
antiparallel to all others in the ab plane but parallel along 
the c-axis), and G-type (every spin is antiparallel to all 
its neighbours). 

From the computations we find that the ground state 
is A-type spin ordered for all values of U. The differ- 
ence in energies between the various magnetically ordered 
structures for LSDA and LSDA+U is reported in TableQ] 
From these total energies it is possible to calculate the ex- 
change constants: we find J\ — 2.6 meV, while J2 = —3.4 
meV for LSDA+U (U = 6.0 and J H = 0.88). These 
quantities can be compared with the exchange constants 
of LaMn0 3 [ll|, where Ji = 9.1 meV and J 2 = -3.1 
meV. 
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FIG. 6: DOS projected on different Cr d-orbitals, calculated 
in LSDA+U with U = 6.0 and J = 0.88 eV for the tetragonal 
structure of KCrF3. Majority/minority contributions to the 
DOS are plotted towards the left/right. The labels ctj, /3, 
label the states 3a; 2 — r 2 , y 2 — z 2 and 3y 2 — r 2 , x 2 — z 2 for the 
different Cr 1 and 2 sites in the ab plane. 




FIG. 7: Contour plot of charge density corresponding to the 
occupied e 9 bands within LSDA+U for U = 6.0 and J = 0.88 
eV for the tetragonal structure of KCrF3 . 



In Fig. [5] the resulting band structure and density 
of states for the tetragonal structure of KCrF3 within 
LSDA+U (U = 6.0 and J = 0.88 eV) are shown. The 
LSDA band gap of 0.485 eV increases to a value of 2.479 
eV (£/=6.0 eV), see Table H 

From the projected density of states we see that within 
LSDA+U around the Fermi level there is a clean distribu- 
tion of the e g states 3x 2 — r 2 /3y 2 — r 2 , depending on the 
Cr-site. A concomitant enhancement of the orbital polar- 
ization is visible in the contour plot of the charge density 
of the occupied e g bands, see Fig. [7] This plot also shows 
that there is an increase in hybrization between the Cr 



e g states and fluoride p states which enhances the total 
band width of the occupied Cr 3d e g bands to about 2.1 
eV. The two-dimensional character of the occupied e g 
bands does not change in the LSDA+U treatment, but 
the dispersion of the empty e g states comes to the fore 
more clearly. From the computations on the cubic phase 
in the next section it will be particularly clear that the 
two-dimensional character of the occupied e g bands that 
is caused by the orbital ordering is also the driving force 
behind the A-type magnetic ordering, as can be expected 
on the basis of the Goodenough-Kanamori [l(| rules for 
superexchange. 
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LSDA+U cubic structure 
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FIG. 8: LSDA+U (17=6.0 and J H =0.88 eV) band structure 
and projected DOS for cubic KCrF3 with Pm — 3m symmetry 
and A-type magnetic ordering. 




FIG. 9: Orbitals in the cubic phase of KCrF3, obtained with 
LSDA+U. Left: orbital ordering of x 2 — y 2 orbitals for C-type 
spin ordering. Right: orbital ordering of 3z 2 — r 2 orbitals for 
A-type spin ordering. Top panels: side view. Bottom panels: 
top view 



HIGH-TEMPERATURE CUBIC PHASE 

In the cubic Pm — 3m structure (a = 4.231783 A) 
the distances between all Cr and neighboring F-ions are 
equal to 2.116 A and the e g states are locally degenerate. 
We find that at the LSDA level, cubic KCrF3 is metallic, 
for the groundstate A-type magnetic structure. Such is 
expected as in the absence of a Jahn- Teller distortion the 
e g band is half filled, even though it is fully spin polarized. 

In LSDA+U a bandgap of 1.72 eV opens up -the band 
structures calculated by LSDA+U (U — 6.0 and J = 
0.88 eV) with the cubic unit cell with A-type magnetic 
ordering is shown in[51 The correlation-induced Mott gap 
is smaller than the charge gap in the tetragonal structure 
because of the absence of the Jahn- Teller distortion. 

Despite the fact that the Jahn- Teller distortions are 
absent in this structure, there is still an orbital order- 
ing which is due to the magnetic exchange, see Fig. [9] 
This exchange-driven orbital ordering can be understood 
in terms of the orbital dependence of the superexchangc 
energy between neighboring Cr-sites. Such a situation is 
described in terms of a Kugel-Khomskii model [23j] . For 
the A-type magnetic ordering we obtain a homogeneous 
orbital occupation of 3z 2 — r 2 states, oriented perpendic- 
ular to the ferromagnetic plane. This is in accordance 
with the Goodenough-Kanamori (lfj| rules for superex- 
change: bonds of occupied 3z 2 — r 2 orbitals on top of 
each other have a large overlap and therefore result in 
antiferromagnetic spin ordering. Within the plane the 
overlap is mainly between occupied 3z 2 — r 2 and empty 
x 2 — y 2 orbitals, causing a ferromagnetic orientation of 
the spins. 

When we consider C-type magnetic ordering the re- 
sulting orbital ordering (see Fig. [9|) is of the homoge- 
neous x 2 — y 2 type -again with antiferromagnetic spin 
orientation for orbitals with lobes pointing towards each 



other (x 2 — y 2 orbitals in the plane) and ferromagnetic 
orientation between orbitals with small overlap. 

The A- and C-type order that we considered are just 
two of many possible magnetic orderings with concomi- 
tant orbital orderings: other configurations can appear 
when the unit cell is doubled or quadrupled in accordance 
with model calculations on the Kugel-Khomskii Hamilto- 
nian These observations of magnetically 

driven orbital ordering (or orbitally driven magnetic or- 
dering, depending on one's point of view), although in- 
teresting from a theoretical perspective, are not directly 
relevant to the experimental situation as at the high tem- 
peratures where the cubic phase is stable no long range 
magnetic ordering is expected. 

By comparing the total energies of the cubic and 
tetragonal phases in the paramagnetic state we can di- 
rectly compute the energy gain in the tetragonal phase 
that is due to the Jahn- Teller distortion alone. We ob- 
tain a value of AEjt = 0.328 eV per unit formula, which 
is comparable to that found in LaMn03 (0.504 eV) 



LOW-TEMPERATURE MONOCLINIC PHASE 

Below 250 K, K&F3 shows a phase transition to a 
monoclinic structure, characterized by a pronounced tilt- 
ing of the CrF@ octahedra. The lattice parameters at 150 
K are a = 5.82642 A, b = 5.83517 A, c = 8.57547 A and 
7 = 93.686° [8]. We performed the total energies calcu- 
lations with a kinetic cutoff energy of 500 eV and used 
the tetrahedron method with Blochl correction, using 90 
irreducible k-points. 

The total energy that we compute in the monoclinic 
phase reveals that the Jahn- Teller distortion is further 
stabilized with an energy gain of 22 meV per Cr with 
respect to the tetragonal phase. Again in the LSDA+U 
calculations the magnetic groundstate is found to be A- 
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LSDA+U monoclinic structure 
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FIG. 10: Monoclinic phase of KCrFs: band structure and 
DOS calculated in LSDA+U with U = 6.0 and J = 0.88 eV. 



type, the Cr moment is 3.85/itt,, the band gap 2.49 eV 
and the orbital ordering is essentially the same as in the 
tetragonal phase (see Fig. fTTj) . Per unit formula the ferro- 
magnetic configuration is higher in energy by 0.0168 eV, 
the C-type configuration by 0.0438 eV and the G-type 
state by 0.0275 eV. 

From this we find slightly larger in-plane and inter- 
plane magnetic couplings of J\ — 2.1 meV and 3 2 = —1.7 
meV, respectively. Similar to LaMnO?; (27j . experimen- 
tally a small spin canting is observed, giving rise to weak 
ferromagnetism in the monoclinic phase [8j. When the 
CuF 6 octahedra are tilted the weak local anisotropy and 
non-local Dzyaloshinskii-Moriya interaction lead to spin 
canting. In a first principles bandstructure calculation, 
this can be taken into account when on top of the present 
computation scheme the relativistic spin-orbit interac- 
tion is included 11 1. 



CONCLUSIONS 

With a set of density functional calculations we de- 
termine the electronic, magnetic and orbital properties 
of KCrF 3 . Our ternary chromium fluoride shows many 
similarities with LaMnOa. From the electronic point of 
view, the bandgap and conduction bandwidth are com- 
parable, although somewhat smaller for the more ionic 
KCrF3 in the orbitally ordered phase. The magnetic 
structure is of the same A-type and the exchange con- 
stants are of the same order of magnitude. The orbital 
ordering in the ferromagnetic planes is identical in the 
two compounds, although the stacking of the ordering 
along the c-axis is different. These properties of KCrF3 
make it a material that is comparable to LaMn03 and 
therefore attractive to investigate for instance orbital ex- 
citations [2 5| . [ 2R [29L |30T | and orbital scattering in pho- 
toemission [3 1|] . Doping the strongly correlated Mott in- 
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FIG. 11: Contour plot of charge density corresponding to the 
occupied e 9 bands below the Fermi level for the monoclinic 
structure of KCrF 3 in LSDA+U. 



sulator KCrF3 with electrons or holes may lead to very 
interesting prospects, as the equivalent manganites show 
an overwhelming wealth in physical properties. If the 
concentration and kinetic energy of the doped carriers 
suffices a melting of orbital ordering is anticipated, es- 
tablishing an orbital liquid phase, changing the electronic 
dimensionality from effectively two to three [i|,[32| ■ In the 
manganites a colossal magneto-resistance is observed in 
the vicinity of such a phase transition. 
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